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We consider the exactly soluble Edwards-Wilkinson Model in one dimension and 
demonstrate explicitly, that it is possible to construct a field 0(r) , that does not 
depend explicitly on time, such that the corresponding time dependent correlation 
function, % q {t) =< <j> q (t)(j>l q (t) > , is dominated at long times by a stretched exponential 

decay. The difference between this and the stretched exponential decay present in 
truly non linear systems is discussed. 



We have argued recently, 1 " 3 that stretched exponential decay is much wider spread 
in condensed matter physics than formerly expected. In fact, we claim, that such a 
form of decay is present in many non linear systems such as, the KPZ dynamical 
systems, ferromagnets at their critical point etc.. These systems have in common the 
property that a disturbance of wave vector q, decays at equilibrium with a 

characteristic decay rate, (O q , that is proportional to q* 1 with ji > 1 . 

The purpose of the present article is to demonstrate explicitly that this is also true 
for linear systems and work out explicitly the form of decay for certain correlation 
function in the Edwards-Wilkinson 4 system. Because the system is exactly soluble 
the full form of the decay function can be written in closed form. 

The Edwards-Wilkinson model is the simplest model describing the landing of 
material and its rearrangement on a growing surface. If we denote the height at a 
point x relative to the spatial averaged height by h(x) , the model is described in one 
dimension by the Langevin equation 

dh 

— (x,t)=vV 2 h(x,t)+T](x,t), (1) 
at 

where the noise, 77 , obeys 

< j](x, t) >= ; < 7](x, t)r](x' ,t')>= 2Dd(x -x')8(t-f) . (2) 

The Langevin equation can be transformed into a Fokker-Planck equation for the 
probability density, P{h } , to find the Fourier transform h , of h(x) . 



1 v 2 

Using the standard transformation, 5 " 8 P = exp[- — ^—j^-h q h_ q ] \j/ , we obtain for y/ a 
Schroedinger like equation 



%■ = -«¥. (4) 
dt 



where the "Hamiltonian" 



H = rD{--^— + l -( V ^fh q h_ q -± V A (5) 
U dhM a 4 D q q 2 D 

H H 



is Hermitian non-negative definite and its only eigenfunction with eigenvalue zero is 
1 v 2 

y/ G ex P[ _ ~52~^"\^- 9 ] • It proves useful to transform to creation and destruction 
operators by writing 



D - 

h q =[— I ]Ha q +al q ) (6) 



and 



3 1 r D ± + 
3^ = -2 [ V ] 

(7) 



where the a' 5 obey the usual Bose commutation relations 



[a,,og = [<,<] = ; [a 9 ,«;] = 5, p . (8) 



In terms of these operators, 

H = Yvq 2 a + q Ct q , (9) 

namely a Hamiltonian of free excitations 

Consider next a general composite field 0(x) . Its time dependent structure factor 

is 

r q (t)=<f_ q m_ q (t)> , (10) 

where the meaning of the average above is that tp q is measured at steady state, the 
system is then allowed to evolve freely and then after a time t , <p_ q is measured. It 
is a standard procedure to show that 



r,(o=<oi0, e->;io> , 



(11) 



here the vacuum 1 > is the normalized state with eigenvalue zero of H . Introducing 
eigenstates of H with momentum q , we obtain the general result, 



r ? (o = £i<oi0 ? i^, J s>i 2 e 



(12) 



where I q, > is a normalized eigenfunction of H that is also an eigenfunction of 
the momentum operator with eigenvalue q and X qp is the corresponding eigenvalue 
of H . 

It is clear from eq. (12), that a necessary condition for a slower than exponential decay 
of T q {t) is that for any q the spectrum A ?/3 is gapless. 

In our case 



where n, is the occupation number of the state £ . The corresponding eigenvalue is 



Consider the state with n q , = n and n e = for any other value of £ . The 



\q,fi >=\{n e }> (such that V £n f = q) , 



(13) 



\ fi ='Yy^n l (such that ^£n e = q) . 



(14) 



corresponding eigenvalue is 



2 




(15) 



Since n can take any integer value it is obvious that the spectrum X q p is indeed 

gapless and the necessary condition for slow decay is fulfilled. 

Now consider for example the case where <j)(x) = h(x) . It is obvious that the 

D _ 2 

decay is not slow but actually T q (t) = — j e vq ' . The reason is that the only matrix 

element of the form < I h q I q , /J > that is not zero is the state with n q = 1 and all 
other n t ' s being zero. If we take (j)(x) = h"(x) (recall that //does not have a zero 
Fourier transform) we arrive again at the conclusion that the decay of T q {t) is 
exponential, because the matrix elements do not couple with states that have a 
vanishingly small "energy". In fact, if (j)(x) = h"(x) the decay will be not slower than 

-v—t 

e " . The conclusion is that in order to obtain a slow decay we must consider such a 
composite field that will couple the vacuum to states with arbitrarily small "energies". 
We consider here the example (j)(x) = exp[iah(x)] , where obviously a is a 

constant that has the dimensions of h~ x . It is straightforward to show that for such a 
composite field 

r 1 9 

<0 9 (O)0_ 9 (O>= \dx exp[--a w(x,t)-iqx], (16) 

2 

where 

t - x 

w(x,t)=<[h(x,t)-h(0,0)] 2 >=D[-Yf( n =) . (17) 

v Vw 



The time dependent correlation function can be written now in the form 



r,(0 = Vw \dy cxp[-a 2 D(-yf(y)] . (18) 



We are interested in the behaviour of T (t) for qV I a 2 D « 1 and Vq 2 t » 1. It is 
easy to show that under these conditions only small values of y are relevant to the 
integration. The conclusion is 



i 



T q {t) = A,\— exp[- 2 qi ], (19) 



where A= f ^-[l-e y ] 

J—oo , 



We see here that the fact that the spectrum is gapless for any q manifests itself in a 
decay that is dominated by a stretched exponential in time, Note, however, that the 
structure is not of the form T q f{(O q t) as is expected in the non linear problem. The 
reason is that in the definition of the field <j)(x) we introduced a dimensional constant 
a and as a result a length, x , into the problem that does not appears in the physics. 
The explicit dependence of the correlation function on the new length x a = v / a 2 D is 
expressed by the following form of T , 



4?/2 

r,(f) = ^-r^i/. 



[^] /2 (ft)„0 /4 exp[- 



(^)(ft) ? /2 ] 



(20) 



This long time decay is very different from the decay of the function < h q (t)h_ q (0) > 

in non linear systems of the KPZ nature or of dynamical (j) 4 theory at the transition , 
that is also found to be dominated by stretched exponential decay. 

Here we have no restriction on the q dependence of the exponential. In the 
non-linear systems, we studied 1 ' 2 ' we found that the non linearity forces that the 
exponential describing the decay of say < h q (t)h q (0) >or analogues of it, to be linear 

in I q I and since the dependence of the exponential on time came only through the 
combination 0) q t , it follows that the stretched exponential factor must be of the form 

exp[-c I I /"^] 9 . 
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